Determination of the Parity of the Neutral Pion via the Four-Electron Decay 
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We present a new determination of the parity of the neutral pion via the double Dalitz decay 
ir° — > e + e~e + e~ . Our sample, which consists of 30 511 candidate decays, was collected from Kl — > 
7r 7r°7r° decays in flight at the KTeV-E799 experiment at Fermi National Accelerator Laboratory. 
We confirm the negative ty° parity, and place a limit on scalar contributions to the ty° — > e + e~e + e~ 
decay amplitude of less than 3.3% assuming CPT conservation. The 7r°7*7* form factor is well 
described by a momentum-dependent model with a slope parameter fit to the final state phase 
space distribution. Additionally, we have measured the branching ratio of this mode to be B(tt° — > 
e+e'e+e-) = (3.26 ±0.18) x 10 -5 . 

PACS numbers: 14.40.Aq, 13.40.Gp 



The parity of the neutral pion has been determined 
indirectly by studying negative pions captured on deu- 
terium [UH]. The observed reactions imply that the w~ 
is a pseudoscalar and that the parities of the tt~ and 
the 7T° are the same. It has long been known that the 
decay 7r° — * 77 in principle offers a direct means of de- 
termining the ir° parity through the polarizations of the 
photons [3J|3]. Given that there are no available meth- 
ods for measuring the polarization of a high-energy pho- 
ton, this measurement has never been performed. How- 
ever, it was soon noted that the double Dalitz decay 
7T° — > e + e - e + e~, which proceeds through an interme- 
diate state with two virtual photons (see Fig. [l~j, is sen- 
sitive to the parity of the pion since the plane of a Dalitz 
pair is correlated with the polarization of the virtual pho- 
ton [51 [5]. This process was studied in a 1962 hydro- 
gen bubble chamber experiment using stopping negative 
pion capture (tt~p — > mr°). That group observed 206 
7T° — * e + e _ e + e _ events and reported that the observed 
distribution of the e + e~ planes was consistent with a 
pseudoscalar pion and disfavored a scalar pion at the 
level of 3.6 standard deviations [7]; this experiment also 
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FIG. 1: Lowest order Feynman diagram for n — » e e~e e~ . 
The direct contribution is shown; a second diagram exists 
with ef and exchanged. 



produced a measurement of the branching ratio of this 
decay, which remains the most precise result to date. 

Using a sample of more than 30 000 tt° — > e + e~e + e~ 
decays, we report new precise measurements of the prop- 
erties of this decay. Our modeling of the decay includes 
for the first time a proper treatment of the exchange 
contribution to the matrix element, and consideration 
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of full 0(a 2 ) radiative corrections. With these advances, 
we have tested for a scalar contribution in the 7r°7*7* 
coupling with a sensitivity of a few percent. We have 
also measured for the first time the momentum depen- 
dence of the form factor in this decay mode. In addition, 
we present a new measurement of the tt° —> e + e~e + e~ 
branching ratio, taking into account radiative effects. 

The most general interaction Lagrangian for the ir° — > 
7*7* transition can be written [5]: 

C oc C^peF^FV® (1) 

where F^ and F pcr are the photon fields, $ is the pion 
field, and the coupling has the form 

C^vpa <xf(xi, x 2 ) [cos (e^vpv 

+ sin^e 1 {,g w g va - g^gup)]- 

The first term in Cp, vpa is the expected pseudoscalar cou- 
pling and the second term introduces a scalar coupling 
with a mixing angle £ and a phase difference S. Nu- 
clear parity violation would introduce a nonzero £, while 
CPT violation would cause the phase S to be nonzero. 
We assume the standard parity-conserving form for the 
7* — > e + e _ conversion. 

The form factor f{x\,X2) is expressed in terms of the 
momentum transfer of each of the virtual photons, or 
equivalcntly the invariant masses of the two Dalitz pairs: 
x\ = ( m e + e - /M„o) 2 ;X2 = (m e + e - /M^o) 2 . In calculating 
the phase space variables for an individual event, there 
is an intrinsic ambiguity in assigning each electron to a 
positron to form a Dalitz pair. Our analysis uses a matrix 
element model that includes the exchange diagrams and 
therefore avoids the need to enforce a pairing choice. 

The 7r°7*7* form factor has been studied previously in 
the decay it — ► e + e~7 [§1[TUJ[TT], where the quantity of 
interest has been the slope parameter a of the first-order 
Taylor expansion f(x,0) = 1 + ax, with x = rrvt + /Ml . 
Here we use a form factor parametrizaton based on the 
model of D'Ambrosio, Isidori, and Portoles (DIP) [T!2] . 
but with an additional constraint that ensures the cou- 
pling vanishes at large momenta |13j . In terms of the 
remaining free parameters, the form factor is: 

, . 1 - fi(l + a)(xi + x 2 ) , , 

fmp{xi,x 2 ;a) = — — — , (3) 

(1 - (ixi)(l - ^x 2 ) 

where \i = M 2 /M 2 w 0.032. In the limit of small x, this 
coincides with the Taylor expansion provided a = —fia. 

The parity properties of the decay can be extracted 
from the angle cf) between the planes of the two Dalitz 
pairs in Fig. [l] where pair 1 is defined as having the 
smaller invariant mass. The distribution of this an- 
gle from the dominant direct contribution has the form 
dT/dcj) ~ 1-Acos(20) + Bsin(20), where A w 0.2cos(2C) 
and B w 0.2sin(2£) cos (5. A pure pseudoscalar coupling, 
therefore, would produce a negative cos(2</>) dependence. 



The 7T° decays used in this analysis are the result of 
fully-reconstructed Kl — > 7r 7r°7r° decays in flight col- 
lected by the KTeV-E799 experiment at Fermilab. The 
E799-II experiment and the KTeV detector are described 
elsewhere [T31 EH]. This analysis relies on two core 
systems of the KTeV detector: a drift chamber-based 
charged particle spectrometer and a cesium iodide (Csl) 
electromagnetic calorimeter. Electrons are identified as 
charged particles whose entire energy is deposited in 
the Csl, while photons are reconstructed from electro- 
magnetic showers in the Csl with no associated charged 
tracks. 

The signal mode, denoted by Kl — > 7r°7r°7rQ D where 
7Tq D refers to 7r° — > e + e~e + e~, has a signature of four 
charged particles identified as electrons and with a com- 
bined invariant mass consistent with the ir° mass, plus 
four photons that are compatible with two additional 
7r°'s. Furthermore, the eight-particle state has an invari- 
ant mass consistent with the Kl and total momentum 
vector in the direction of the kaon line of flight. 

The branching ratio measurement, which we describe 
here first, makes use of a normalization mode in which 
two pions decay via ir° — > e + e~7 and the third 7r° — > 
77. This "double single-Dalitz" mode, denoted Kl — > 
7r°7Tj! ) 7rQ where 7Tq refers to 7r° — > e + e~7, has the same 
final state particles as the signal mode and is again iden- 
tified by finding the proper combinations of particles to 
make three pions with a total momentum consistent with 
the kaon. The similarity of these modes allows cancel- 
lation of most detector-related systematic effects in the 
branching ratio measurement, but also allows each mode 
to be a background to the other. 

Radiative corrections complicate the definition of the 
Dalitz decays in general. We define the signal mode 
7T° — > e + e~e + e~ to be inclusive of radiative final states 
where the squared ratio of the invariant mass of the four 
electrons to the neutral pion mass x 4e = (M 4e /M w o) 2 
is greater than 0.9, while events with x 4e < 0.9 (ap- 
proximately 6% of the total rate) are treated as 7r° — > 
e + e~e + e~7- For normalization, the decay ir° — > e + e _ 7 
is understood to include all radiative final states, for con- 
sistency with previous measurements of this decay |16j . 
Radiative corrections in this analysis are taken from an 
analytic calculation to order 0(a 2 ) |H]. 

Other final states of the Kl — > ir a ir°ir° decay can 
become backgrounds to either the signal or normaliza- 
tion mode if one or more photons convert to an e + e~ 
pair in the detector material: Kl — > 7r 7r 7TQ where one 
of the five photons converts, or Kl — > tt°tt°tt — > 67 
where two photons convert. These modes again have 
the same final state as the signal, but can be distin- 
guished statistically since the externally produced pairs 
tend to have smaller invariant masses than those from 
internal conversions. The most significant of these back- 
grounds is Kl —> 7r 7r 7Tp with one external conversion 
in material. The photon must convert upstream of the 
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FIG. 2: Invariant e + e~e + e~ mass for data events identified 
by the preferred pairing as Kl — * 7r°7r°7rfj D (upper panel) and 
Kl — > 7r 7rp7rp (lower panel). 



first drift chamber for the resulting tracks to be recon- 
structed. The material in this region sums to 2.8 x 1CP 3 
radiation lengths. With five photons available, the prob- 
ability of one converting is 1.08%, close to the single- 
Dalitz branching ratio. The distinguishing characteristic 
of these events is the small value of the e + e~ invariant 
mass, or similarly, the small value of the opening angle 
of the pair. Requiring a track separation at the first drift 
chamber of greater than 2 mm removes 99.74% of the 
remaining simulated background while preserving 74.3 % 
of signal and 72.7% of normalization events. 

The final selection criterion separates Kl 7r ° 7r ° 7r DD 
from Kl — > 7r ° 7r D 7, D events. This is accomplished by a \ 2 
formed of the three reconstructed tt° masses. This serves 
to identify the best pairing of particles for a given decay 
hypothesis, as well as to select the more likely hypothesis 
of the two. The event is tagged as the mode with the 
smaller x 2 , which is further required to be less than 12 
(with three degrees of freedom) . This technique correctly 
identifies more than 99.5% of events (Fig. [2]). 

The final event sample contains 30 511 signal candi- 
dates with 0.6% residual background and 141 251 nor- 
malization mode candidates with 0.5% background (de- 
termined from the Monte Carlo simulation). The back- 
ground in the signal sample is dominated by mistagged 
events from the normalization mode. 



The branching ratio is measured from the ratio of re- 
constructed signal mode events to normalization mode 
events. This ratio must be corrected by the ratio of ac- 
ceptances, which has been determined using a detailed 
Monte Carlo simulation of the beam distribution and de- 
tector response. The resulting double ratio is directly re- 
lated to the branching ratio B eeee of the 7r° — > e + e~e + e~ 
mode normalized to the square of the branching ratio 



Beeee ' B^^ 

B 2 



N(K L 



'DDJ 



(K L 



'd"d; 



iV(X L ^7r%^) e(K L 



r u 7if )n ) ' 



. (4) 

where N is the number of events and e is the combined 
geometric acceptance and detection efficiency for a given 
mode. 

The statistical error on the ratio in Eq. [4] is 0.62%. 
Systematic errors on the efficiencies were determined 
through data studies as well as variations in the param- 
eters of the Monte Carlo simulation. Because the final 
state particles in the signal and normalization mode are 
the same, detector-related quantities substantially cancel 
in the ratio, which is generally insensitive to the details 
of the simulation. The dominant systematic errors came 
from variation of the analysis cuts (0.21%) and Monte 
Carlo simulation statistics (0.25%). Other systematic er- 
rors were from uncertainties in the amount of material 
in the spectrometer (0.15%), uncertainty in the back- 
ground levels in the two samples (0.15%), modeling of the 
drift chamber resolutions (0.11%), and radiative correc- 
tions (0.04%). The total systematic error on the relative 
branching ratio is 0.41%. 

The final result for the ratio of decay rates is: 
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=0.2245 ± 0.0014(stat) ± 0.0009(syst). 



(5) 



The 7T° — > e + e e + e branching ratio can be calculated 
from the double ratio using the known values -B 77 = 
0.9880 ± 0.0003 and B eei = (1.198 ± 0.032) x 10" 2 [17]. 
This yields B*>° e 9 = (3.26 ± 0.18) x 10~ 5 , where the er- 
ror is dominated by the uncertainty in the ir° — > e + e _ 7 
branching ratio. Using our radiative corrections model 
[5] to extrapolate to all radiative final states, we find: 
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7 - =0.2383 ± 0.0015(stat) ± 0.0010(syst), 

(6) 

(3.46 ±0.19) x 10" 5 . Our branching ratio 
result is in agreement with previous measurements [7]- 

The parameters of the 7r°7*7* coupling are found by 
maximizing an unbinned likelihood function composed 
of the differential decay rate in terms of ten phase-space 
variables. The first five are (xi,X2,yi,y2,^>), where x±, 
X2, and cj> are described above and the remaining vari- 
ables yi and yi describe the energy asymmetry between 
the electrons in each Dalitz pair in the 7r° center of mass 
[8]. The remaining five are the same variables, but cal- 
culated with the opposite choice of e + e~ pairings. The 
likelihood is calculated from the full matrix element in- 
cluding the exchange diagrams and 0{a 2 ) radiative cor- 
rections. 
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The fit yields the DIP a parameter and the (complex) 
ratio of the scalar to the pseudoscalar coupling. For rea- 
sons of fit performance, the parity properties are fit to the 
equivalent parameters k and 77, where k + irj = tan£e 1<5 . 
The shape of the minimum of the likelihood function in- 
dicates that the three parameters a, k, and 77 are uncor- 
rected. Acceptance-dependent effects are included as a 
normalization factor calculated from Monte Carlo simu- 
lations. 

Systematic error sources on a and k are similar to those 
for the branching ratio measurement. The dominant er- 
ror is due to variation of cuts, resulting in a total system- 
atic error of 0.9 and 0.011 on a and k respectively. For 
the 77 parameter, the primary uncertainty results from 
the resolution on the angle <fi between the two lepton 
pairs, which produces an effective flattening of the angu- 
lar distribution without inducing a phase shift. The fitter 
interprets this as a small scalar contribution with a phase 
difference of 90 degrees, and therefore a larger value of 
77, particularly for 77 « 0. This behavior was studied with 
Monte Carlo simulation and a correction was calculated. 
The uncertainty on this correction results in a systematic 
error of 0.031. 

The distributions of x\ and 22, overlaid with the Monte 
Carlo simulation, are shown in Fig. [3] The cf> distribution 
is shown in Fig. [4] For plotting the data a unique pair- 
ing of the four electrons is chosen such that X\ < X2 and 
the product x%X2 is minimized: this choice represents the 
dominant contribution to the matrix element. It is clear 
that the pseudoscalar coupling dominates, as expected, 
with no evidence for a scalar component. The distribu- 
tions of all five phase space variables agree well with the 
Monte Carlo simulation. 

The final results for the three parameters are a = 
1.3 ± l.O(stat) ± 0.9(syst), k = -0.011 ± 0.009(stat) ± 
O.Oll(syst), and 77 = 0.051 ± 0.026(stat) ± 0.031(syst). 
The DIP a parameter is related to the standard slope 
parameter by a = -0.032a, yielding a = -0.040 ± 0.040. 
This result is in agreement with recent direct measure- 
ments. 

The parameters k and r\ are transformed into limits 
on the pseudoscalar-scalar mixing angle C under two hy- 
potheses. If CPT violation is allowed, then the limit is 
set by the uncertainties in 77, resulting in C, < 6.9° at 
the 90% confidence level. If instead, CPT conservation 
is enforced, 77 must be zero, and the limit derives from 
the uncertainties on k, resulting in £ < 1.9°, at the same 
confidence level. These limits on £ limit the magnitude 
of the scalar component of the decay amplitude, relative 
to the pseudoscalar component, to less than 12.1% in the 
presence of CPT violation, and less than 3.3% if CPT is 
assumed conserved. The limits on scalar contributions 
apply to all ir° decays with two-photon intermediate or 
final states. 

This analysis confirms the negative parity of the neu- 
tral pion with much higher statistical significance than 
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FIG. 3: Distribution of the kinematic variables X\ and X2 
for signal event candidates (points) and signal Monte Carlo 
simulation with best fit form factor parameters (histogram). 




FIG. 4: Distribution of the angle (f>, in units of n, between 
the planes of the two e + e~ pairs. The solid histogram shows 
the Monte Carlo expectation for negative parity. 



the previous result, and places tight limits on nonstan- 
dard scalar and CP T- violating contributions to the ir° — > 
e + e~e + e - decay. We have also measured the momentum 
dependent form factor in this decay for the first time, and 
made the first improvement in its branching ratio since 
1962. This measurement is limited at present by the cur- 
rent large uncertainty in the branching ratio of the single 
Dalitz decays used for normalization, but we expect that 
uncertainty to be reduced in the near future at which 
point the present measurement can be recalculated using 
the more precise double ratio measurement. 
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